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INTRODUCTION 


It  can  be  said  that  a majority  of  modern  high  thrust  propulsion 
nozzles  obtain  their  thrust  by  means  of  a reactive  force  to  the  ex- 
pelled fluid.  Prediction  and  maximization  of  the  performance  of  these 
devices  require  that  the  flow  field  be  accurately  analyzed  in  the 
three  flow  stages,  naniely  the  subsonic,  the  transonic,  and  the  super- 
sonic regimes. 

In  the  present  investigation,  the  transonic  flow  region  of  a pro- 
pulsion nozzle  is  studied.  In  this  region,  high  pressure  fluid  is 
received  as  a result  of  some  combustion  process.  The  second  stage 
receives  potential  energy,  associated  with  static  fluid  temperature, 
and  relinquishes  kinetic  energy,  associated  with  translational  fluid 
velocity.  Thus,  in  the  transonic  region,  an  increase  in  linear 
momentum  is  realized.  The  outflow  of  this  stage  is  supersonic  and, 
in  fact,  provides  initial  conditions  for  the  third  stage  along  which 
the  flow  is  everywhere  supersonic. 

To  this  date  there  has  been  abundant  literature  published  on  the 
subject  of  transonic  nozzle  flow.  One  general  method  of  solution, 
and  the  method  of  solution  presented  herein,  is  a series  expansion 
technique.  Generally,  the  expansion  approach  involves  a parametric 
solution  to  a particular  set  of  governing  flow  equations  where  the 
expansion  parameter  can  be  directly  related  to  the  flow  geometry. 

One  of  the  pioneers  of  the  expansion  method  of  solution  was 
Sauer.  In  his  original  paper  (1)  the  governing  equations  are  reduced 
to  a convenient  form  by  means  of  a perturbation.  The  flow  velocity  is 
then  expanded  in  a power  series  as  a function  of  the  radial  coordinate. 
The  coefficients  of  the  series  are  obtained  by  substituting  the  series 
into  the  perturbed  governing  equations.  These  coefficients  turn  out 
to  be  functions  of  the  centerline  velocity. 

A different  approach  to  the  expansion  technique  was  used  by 
Hal  1.(2)  In  this  report,  the  velocity  components  along  the  region  of 

^^^Sauer,  R.,  "General  Characteristics  of  the  Flow  Through  Nozzles  at 
Near  Critical  Speeds,"  NACA  TM  No.  1147  (1947). 

(2) 

'Hall,  I.M.,  "Transonic  Flow  in  Two-Dimensional  and  Axial ly-Symnetric 
Nozzles,"  Quart.  Journ.  Mech.  and  Applied  Math.,  Vol . XV,  Pt.  4 (1962) 
pp.  487-508. 
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the  throat  are  calculated  by  means  of  an  expansion  in  inverse  powers 
of  R,  the  radius  of  curvature  of  the  nozzle  profile  at  the  throat.  It 
can  be  shown(3)  that  Sauer's  solution  is  a special  case  of  Hall's 
solution  when  terms  of  order  R-2  are  neglected.  Two  innovative  aspects 
to  Hall's  work  made  it  a unique  contribution  to  the  literature.  One 
concept  is  to  expand  the  wall  contour  in  terms  of  a McLaurin  series 
and  the  other  is  to  employ  a "coordinate  stretch".  Both  of  these  ideas 
are  used  in  the  following  investigation. 

It  is  observed  that  the  nozzle  boundary  is  not  a constant  coordi- 
nate line  in  cylindrical  coordinates  and  therefore  must  be  evaluated 
as  a power  series.  So,  the  wall  boundary  condition  is  never  exactly 
satisfied  in  cylindrical  coordinates.  Also,  the  radial  velocity 
component  is  proportional  to  the  boundary  slope,  which  can  become 
quite  large  in  the  throat  region  of  nozzles  having  small  radii  of  cur- 
vature. Hence,  Kliegel  and  Levine  (4)  speculated  that  the  accuracy 
of  the  solution  could  be  improved  by  considering  a solution  in  a 
coordinate  system  where  the  flow  axis  and  wall  are  both  coordinate 
lines.  This  suggests  a toroidal  coordinate  system  for  nozzles  having 
throats  whose  cross-section  is  circular.  It  happens  that  in  toroidal 
coordinates,  the  solution  now  appears  in  powers  of  1/(R  +1).  In 
this  system,  the  convergent  properties  of  the  expansion  solution  are 
in  fact  superior.  In  addition,  in  the  limit  of  large  radii  of  curva- 
ture, the  flow  field  of  Kliegel  and  Levine  reduces  to  Hall's,  as  it 
must.  Finally,  it  is  noted  C^)  that  comparisons  between  theory  and 
experiment  are  excellent. 

Another  expansion  method  which  uses  an  inverse  approach  for 
solving  the  transonic  flow  problem,  was  proposed  by  Hopkins  and  Hill.^  ' 
In  their  development,  a transformation  to  a velocity-potential,  stream- 
function  coordinate  system  was  used  in  simplifying  the  governing 
equations.  Four  partial  differential  equations  result  from  the  basic 
flow  characterization.  The  partial  differential  equations  are  solved 
by  means  of  writing  the  unknowns  in  the  form  of  a power  series;  the 
coefficients  being  functions  of  velocity  potential  only.  The  two 
authors  make  use  of  the  solution  to  analyze  annular  type,  unconvention- 
al nozzles. 


u) 

'May,  R.J.  Jr.,  Thompson,  H.D.,  and  Hoffman,  J.D.,  "Comparison  of 
Transonic  Flow  Solutions  in  C-D  Nozzles,"  Air  Force  Aero  Propulsion 
Laboratory  Report  No.  AFAPL-TR-74-110,  (October  1974). 

(4) 

' 'Kliegel,  J.R.,  and  Levine,  J.N.,  "Transonic  Flow  in  Small  Throat 
Radius  of  Curvature  Nozzles,"  AIAA  Journal,  Vol.  7,  No.  7 (July 
1968)  pp.  1375-1378. 

(5) 

' 'Hopkins,  D.F.,  and  Hill,  D.E.,  "Transonic  Flow  in  Unconventional 
Nozzles,"  AIAA  Journal,  Vol.  6,  No.  5 (November  1967),  pp.  838-842. 
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Until  now,  almost  all  of  the  transonic  flow  analyses  have  dealt 
with  either  an  axisynmetric  or  planar  nozzle  geometry.  However,  in 
present  design,  it  is  sometimes  necessary  to  use  non-circular  pro- 
pulsion nozzle  to  facilitate  the  nozzle,  airframe  integration.  Ivanov 
and  Ryl'ko(6)  consider  a three-dimensional  transonic  flow  problem  for 
an  elliptic  geometry  by  writi’q  the  system  of  governing  equations  (gas 
dynamic  equations)  in  integra.  form.  The  equations  are  solved  numeri- 
cally! 7, 8)  using  a difference  scheme  of  continuous  calculations  of 
the  discontinuities  of  the  solutions  to  the  gas  dynamic  equations.  To 
accomplish  this,  one  fourth  of  the  flow  field  (from  symmetry  consider- 
ations) was  divided  up  into  3200  elemental  cells.  The  initial  distri- 
bution of  flow  variables  in  the  nozzle  was  taken  from  the  one-dimen- 
sional computations.  Various  ratios  of  semi-major  to  semi-minor  axes 
were  chosen  and  numerical  results  for  the  flow  distribution  are  given. 

This  report  concerns  the  analytical  solution  of  the  transonic  flow 
in  a converging-diverging  nozzle  having  an  elliptical  cross-section 
and  follows  closely  the  M.S.  Thesis  of  Epstein  (9).  A power  series 
expansion  technique  similar  to  that  of  Hall (2)  is  employed.  The 
procedure  consists  of: 

1.  Writing  the  governing  equations  and  boundary  conditions  in 
a suitable  non-dimensional  form. 

2.  Writing  the  velocity  components  as  a power  series  in  terms  of 
an  expansion  parameter  e. 

3.  Determining  the  first  order  governing  equations  and  boundary 
conditions  by  substituting  the  power  series  form  for  the 
velocity  components  into  the  governing  equations. 

4.  Assuming  an  appropriate  (consistant)  form  for  the  solution 
to  the  first  order  governing  equations. 


^^^Ivanov,  M.  Ya.,  and  Ryl'ko,  O.A.,  "Analysis  of  Transonic  Flow  in 
Elliptical  Nozzles,"  IZV.  AKAO.  NAUK  SSSR,  MEKH.  ZHID.  GAZA,  No.  3 
(September  1972),  pp.  161-163. 

^^^Godunov,  S.K.,  "Difference  Method  of  Numerical  Solution  of  Differ- 
ent Problems  of  Fluid  Mechanics."  Matem.  Sb.,  Vol . 47(89),  No.  3 
(1959),  pp.  271-306. 

^^^Godunov,  S.K.,  Zabrodin,  A.V.,  and  Prokopov,  G.P.,  "Difference 
Scheme  for  Two-Dimensional  Unsteady  Problems  of  Gas  Dynamics  and  the 
Calculation  of  Flows  with  Detached  Shock  Waves,"  Zh.  Vychislit.  Matem. 
i Matem.  Fiz.,  Vol.  1,  No.  6 (1961),  pp.  1020-1050. 

^^^Epstein,  M.A.,  "Transonic  Flow  in  Nozzles  Having  an  Elliptical 
Cross  Section,"  M.S.  Thesis,  Purdue  University,  December  1976. 
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5.  Using  tne  boundary  conditions  to  solve  for  the  various 
constants  that  appear  in  the  solution. 

In  the  first  appendix  of  the  work  done  by  May,  Thompson,  and 
Hoffman, (3)  a general  transonic  solution  for  an  axisyinnetric  nozzle 
is  presented.  The  relation  between  the  work  conducted  herein  and 
that  of  reference  (3)  is  discussed  in  Section  V. 
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SECTION  II 


GOVERNING  EQUATIONS 


4. 


The  nozzle  flew  's  considered  to  be  irrotational , inviscid,  and 
steady.  Orthogonal  cylindrical  coordinates  (x,r,0)  are  used,  and  the 
origin  of  the  coordinate  system  is  taken  at  the  center  of  the  throat. 
X is  measured  in  the  downstream  direction  with  r measured  radially 
and  0 measured  tangentially. 


The  governing  equations  are  the  Euler  equations, 

-VP  = pDV/Dt  (1) 

and  the  continuity  equation, 


^ + pV  • V = 0 


(2) 


~y  ~y 

where  V = (ir-Vp  + ievg  + ix^x)’  Vx»  and  vq  are  the  velocity  compo- 
nents in  the  x,  r,  and  0 directions  respectively;  D/Dt  is  the  substan- 
tial derivative  given  by 


Dt 


3r 


r 30 


+ V ^ 
X oX 


(3) 


for  steady  flow,  and  pis  the  density.  Also,  the  "del"  operator  in 
cylindrical  coordinates  is  written: 


7 = 1 T T ^ 

'r  3r  'o  r 30  x 3x 


(4) 


where  ir,  ie  and  i^^  are  unit  vectors  in  the  coordinate  direction 
Now  the  definition  of  the  speed  of  sound,  namely, 

2 


dp 


= a 


(5) 


can  be  written  for  an  irrotational,  steady,  inviscid  flow  in  the  form. 


DP  = a2 
Dt  Dt 


(6) 


Using  equation  (6)  to  eliminate  density  derivatives  from  equation  (2) 
and  expanding  the  result  yields 
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3r  r r r 30  3x 
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r rir 


— ^^-  + V ^ = 0 
r 30  X Dx 


(7) 


where  a is  the  sound  speed  and  P is  the  pressure. 

The  next  step  is  to  use  the  Euler  equations  to  eliminate  the 
pressure  from  equation  (7).  In  cylindrical  coordinates  the  Euler 
equations  can  be  written 


3r 

i = 

r 30 

3P 


- — - elv  -jT  ■*  TW*\W 


9v 


'^r  3r  r 30  '^x  3x 
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9x 


= P 


3v 


''r  3r 


v„  3v 

X . 0 X 


3v 


+ V — - 

r 30  X 3x 


(8) 

(9) 

(10) 


Equations  (8)  through  (10)  are  substituted  into  equation  (7)  so  that 
the  pressure  terms  are  eliminated: 


3V 

^ \l  \l 

[ 1 ^ 

3r 

[ rj  r 30 

r "oj 

3x 

''xJ 

''r''e| 

[ r 30  3r  j 

pv 

- V^x  [3x 

3v 

+ 

3r 

- v,,v 

0 X 

i!!i 

i^r  30 

^''ol 

r r 

Now  for  a steady,  irrotational , inviscid  flow  the  Euler  equations  can 
be  replaced  by  the  irrotational ity  condition,  namely 

V X V = 0 (12) 

Now,  expanding  the  irrotational ity  condition  in  cylindrical 
coordinates  yields 


30 

3x 

(13) 

■3Vr 

3x 

^"x 

W 

(14) 

"''r. 

3v 

3(rVQ) 

(15) 

30 

^ 3r  6 

3r 
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Using  equations  (13)  through  (15)  in  equation  (11)  yields: 


8v 
r 

8r 


fa^-v^l 
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9v 
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r 36 
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9x  '^x. 

V v' 
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3v 
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r X dr 


- 2v„v 


0 X 9x 


= 0 


(16) 


The  nozzle  geometry  near  the  throat  is  assumed  to  be  elliptical 
in  cross-section  so  that  the  boundary  geometry  is  described  by  the 
equati on 


2 2 2-2 
r cos  9 ^ r sin  9 

a^(x)  B^(x) 


(17) 


where  a and  3 are  even  functions  of  x such  that  the  nozzle  is  synmetric 
about  X = 0 in  the  throat  region.  Writing  a and  3 in  a power  series 
for  smal 1 val ues  of  x. 


a(x)  = Aq  + A^x^  + 0(x'^) 

(18) 

3(x)  = Bq  + B2X^  + O(x^) 

(19) 

where  0(x”^)  denotes  terms  of  order  and  higher. 


It  is  convenient  to  nondimensionalize  the  problem  at  this  point. 
The  length  Aq  is  used  to  nondimensionalize  all  distances  and  the 
critical  sound  speed,  a*,  is  used  to  nondimensionalize  velocities. 
Thus , 


r = and  X = ^ 
0 0 


(20) 


0 is  already  nondimensional  since  it  is  measured  in  radians.  Similarly, 


V = 

''r  - a*  ’ 


- _ 0 
''0  = a^ 


and  v^  E ^ 


(21) 


For  the  adiabatic  flow  of  an  ideal  gas  Bernoulli's  equation  can 
be  written 


,2  . rS  ..2  . ri 


2 2 2 

V ''0  ''x 


(22) 


Substituting  equation  (22)  as  well  as  the  nondimensional  coordinate 
and  velocity  components  into  equation  (16)  and  simplifying  yields: 
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= 0 


Using  the  nondimensional  expressions  for  the  velocity  components 
and  coordinates  in  equations  (13),  (14),  and  (15)  gives: 
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(24) 
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3v„ 

r 

3x 

X 

'dr 

(25) 

3v„ 

9v 

r 

30  ’ “ 

r — + 

9r  e 

(26) 

Equations  (23)  through  (26)  are  the  nondimensional  governing 
equations.  The  procedure  now  is  to  write  the  dimensionless  velocity 
components  in  terms  of  an  expansion  parameter,  e,  where  e is  small 
compared  to  one.  Then,  since  consideration  is  to  be  limited  to  a 
region  close  to  the  throat,  v^  and  vq  are  small  and  is  approximately 
equal  to  unity.  The  assumed  perturbation  of  the  coordinate  velocities 
are: 


= 1 " "%1 

* * ■■■ 

(27) 

3/2- 

- ^ Vl 

(28) 

3/2- 

- " ''01 

(29) 

These  particular  powers  of  c are  used  so  that  the  boundary  conditions 
can  be  satisfied  and  so  that  a reduction  to  the  special  case  of 
axisymmetric  flow  is  possible. 

In  order  that  the  irrotationality  equations  (24),  (25),  and  (26) 
be  satisfied,  the  x-coordinate  must  be  transformed  to 

X = (30) 
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PT 


Substituting  equations  (27)  through  (30)  into  the  governing 
equations  (equations  (23)  through  (26))  and  retaining  only  the  lowest 
order  powers  of  epsilon  yields 


3r 

1 

r 

30 

^''zl  1 

(31) 

N 

> 

^''el 

(32) 

30 

r 

3z 

'Vi . 

3v 

zi 

(33) 

3Z 

„ ^''oi 

"Vi 

(34) 

3r 

+ 

< 

O 

II 

30 

For  convenience,  the  bars  over  the  velocities  in  the  three  coordinate 
directions  and  over  r and  z have  been  dropped.  Equations  (31)  through 
(34)  are  the  first  order  governing  equations  for  the  present  problem. 


i 

I 

1 

i 


I 

I 
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SECTION  III 


BOUNDARY  CONDITIONS 


! I 


The  first  boundary  condition  arises  from  the  fact  that  the 
velocity  component  normal  to  the  nozzle  boundary  is  zero.  That  is. 


V • n = 0 

where  n is  a unit  normal  to  the  nozzle  surface. 


(35) 


To  write  this  condition  in  a suitable  form  for  the  present 
problem,  it  is  necessary  to  recall  the  expression  for  the  nozzle's 
geometry,  equation  (17).  This  can  be  nondimensionalized  as  follows. 
Let 


a(x)/A^  = a(x)  h 1 + A^A2X 


-2 


(36) 


and 


Then, 


B 2 

3(x)/A^  5 3(J)  h ^ + A^B^i^ 

2 2 2.2 
r cos  9 r sin  9 _ , 

a^(x)  6^(x) 


(37) 


(38) 


^2  3nd  B2  can  be  real  ted  to  the  expansion  parameter  e_and  will  be  dis- 
cussed later.  Again,  for  convenience,  the  bars  over  a(x)and  B(x)  will 
be  dropped. 

For  a wall  contour  defined  by  r = h(x,0),  a McLaurin  expansion 
about  X = 0,  holding  for  all  theta  yields 

r = h(x,e)  = h(O,0)  + m'x  + ^m"x^  + ^m'"x^+  ...  (39) 


where 


. _ 3b| 


m = 


ax  x=o* 


- Ol 
” 2 
ax 


x=0 


33. 

.Ill  _ h 

~ 3 

3x 


x=0’ 


(40) 


10 


From  equation  (38)  and  from  the  definition  of  the  wall  contour. 


r = h(x,0)  = a(i/[B^cos^0  + uSin^O 


1/2 


(41) 


Then, 


and 


m 


9^h 

3x2 


. _ 9h| 

O 

II 

o 

II 

X 

■ 3x 

1 2A„ fAoB^cos20  + B.A^sin^ol 

1 _ 0 1 Z 0 2 0 1 

[b^cos^o  + A2sin20 
1 0 0 ; 

3/2 

(42) 


(42) 


The  expansion  of  higher  order  derivatives,  i.e.,  m' ' ' , m^'^,  etc.  re- 
sults in  a series  with  terms  of  higher  order  in  epsilon  and  are  there- 
fore neglected. 


Now,  returning  to  equation  (38),  observe  that  it  can  be  written 


as 


? 9 7 9 9 9 9 9 

B r COS  0 + a'^r  Sin  0 - a 8 = 0 = f(r,0,x) 


(44) 


where  f is  a function  of  r,  6,  and  x and  represents  the  nozzle  surface. 
But, 


Vf  = cn 


(45) 


where  c is  a constant.  In  addition. 


3r  V r 30  ^0  3x  x 


(46) 


Thus,  equation  (35)  may  be  written  as; 


V 31+  lilf  + V 0 

3r  r 30  ''x  3x  ^ 


(47) 


Next,  the  quantities  v , v.,  and  v„  will  be  expanded  about 
r = h(O,0)  H h„;  ® ^ 


= v^(x,h,9)  = v^(x,h^,0)  + (h-h^) 


!!r 

3r 


(h-h^)2  3^v^ 
3 r~^ 


X , hp , 0 


(48) 


x,h^,0 
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V„  = v,,(x.h,.')  - Vy(x.h^.o)  + (h-h^)  ^ 


x.h^.e 


2: 


i)f3^ 


X » j 0 


(49) 


= v^(x,h,y)  = v^(x,h^,0)  + (h-h^)  y-^- 


(h-h  3^v 
0'  X 


x,hg,0 


2: 


ar 


x,h^,0 


(50) 


Then , 


Now,  recall  that  the  wall  contour  was  defined  by  r = h(x,0), 


r - h(x,0)  = 0 = f(r,0,x) 

Differentiating  this  equality  gives  the  following: 

31-1  3f^  ^ ;.na3f_  3h 

3r  ’ 30  ■ 30  ’ 3x  ' 3x 


(51) 


(52) 


This  factor  1/r  can  be  ezpanded  in  terms  of  the  boundary  by  means  of 
an  expansion  along  the  nozzle  wall.  Observe  that 


1 _ 1 


r ■ h^-^6  ■ hJUd/hJ  - h^ 


1-  Jl-  + ^ 
. 2 

■ ° ^ 


(53) 


Then,  equation  (47)  becomes 


V - V 
r ''0  30  h. 


1 - ^ 
h . 2 


- V 0 

X 3x 


(54) 


The  expansions  for  v^,  Vy,  and  are  substituted  from  equations  (27), 
(28)  and  (29)  and  higher  order  terms  are  dropped.  Thus, 


- • 0 (55) 
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But  can  be  found  directly  from  equation  (39): 

oX 


= m'  + m' ‘x  = m' 
9x 


(56) 


since  m'  was  found  to  be  zero.  Therefore,  equation  (55)  becomes 


c^^^v^,(2.h„,e)  - c3/2,^,(z.h„.6)  II  C 


- c^'^^m"v^^(z,h^,0)  z = 0 


(57) 


where  m' ' is  given  by  equation  (43). 

In  order  that  the  terms  of  equation  (57)  be  of  the  same  order  in 
e and  reduce  to  the  axi symmetric  case  the  following  relationship 
between  and  c is  assumed. 


= e/2;  = ke/2. 


(58) 


where  k is  a constant.  The  e/2  terms  are  essential  for  a reduction 
to  the  axi symmetric  case  but  the  selection  of  constant  k is  arbitrary. 
Moreover,  it  is  observed  that  k has  the  effect  of  restricting  the 
nozzle  geometry  by  fixing  the  ratio  between  nozzle  semi -major  and  semi- 
minor  axes. 

By  the  relationships  of  equations  (58)  m"  is  of  order  e and 
equation  (57)  reduces  to 


.3/2 


Vl 


(z.hQ,e) 


'ol 


(z.h^,©) 


9h 

96 


m'  'z 

e 


= 0 


(59) 


which  is  the  first  boundary  condition. 

A second  boundary  condition  arises  from  the  fact  that  the  radial 
velocity  at  the  nozzle  centerline  ( r = 0)  must  be  zero.  That  is, 

v^^(z,O,0)  = 0 (60) 


Finally,  consideration  of  the  tangential  velocity  component  offers 
two  additional  boundary  conditions.  Since  the  two  geometrical  planes 
of  syrmetry  are  chosen  as  the  planes  of  symmetry  for  the  tangential 
velocity  component,  this  velocity  component  vanishes  at  0 = 0,  tt/2, 

IT,  and  3tt/2  for  all  values  of  r and  z.  Due  to  elliptical  symmetry  only 
one-fourth  of  the  flow  region  needs  to  be  considered.  So, 
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(61) 


VQ,(z,r,0)  = 0 

(z,r,7i/2)  = 0 (62) 

Equations  (59),  (60),  (61),  and  (62)  are  the  four  first  order 
boundary  conditions. 


1 
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SECTION  IV 


PROBLEM  SOLUTION 


From  Section  II  the  first  order  governing  equations  are: 

3v  , 9v  9v  V 

^ ^ r 9^  ■ v^(y+1  ) + — - 0 


(31) 


9V. 

90 


9Z 


(32) 


9v^  9v^ 

IF  " ^ 


(33) 


±_  (rv  ) = — - 
9r  ^ q'  96 


(34) 


and  from  Section  III  the  first  order  boundary  conditions  are: 


v^(z,hQ,e) 


VQ(z,ho,0) 


(59) 


v^(z,O,0)  = 0 
Vg(z,r,0)  = 0 
Vg(z,r,Tr/2)  = 0 


(60) 

(61) 

(62) 


where  the  subscript  one  has  deen  dropped  from  all  velocities  since  con- 
sideration is  limited  only  to  first  order.  Note  that  the  governing 
equations  and  boundary  conditions  for  axi symmetric  flow  are  a special 
case  of  the  present  formulation  and  are  obtained  by  setting  V0=9/96=O. 
Further,  it  will  be  shown  in  Section  V that  the  solution  to  the 
present  three-dimensional  transonic  flow  problem  reduces  to  the  axisym- 
metric  solution. 
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To  begin  the  solution,  a series  of  the  form 

dE^(r.O) 

+ E^(r,0)  + z (63) 

is  assumed  for  the  radial  velocity  component  where  C-|  is  a constant  and 
El  and  E2  are  functions  of  r and  0.  is  chosen  in  this  fashion  so 
that  it  can  be  represented  as  a linear  yet  completely  general  function 
of  r,  0,  and  z,  and  at  the  same  time  reduce  to  the  axisymmetric  case 
when  the  0 dependence  is  surpressed. 

Using  this  representation  of  v^,  it  is  seen  that  in  order  to 
satisfy  equation  (33)  the  axial  component  of  velocity,  Vp^,  must  be  of 
the  form 


v^  = E2(r,0)  + D^(0)  + D2(0)z 


(64) 


where  D-j  and  D2  are  functions  only  of  0. 

Similarly,  in  order  that  equation  (32)  be  satisfied. 


V 


0 


£ 

r 


' 3E2(r,0)  3D,(0)' 

30  ^ 30 


1 302(0) 

f 30 


2 + E3(r,0) 


(65) 


r 


where  E^  is  a function  of  r and  6. 

Now,  when  equations  (63),  (64),  and  (65)  are  substituted  into 
equation  (31),  one  obtain  a quadratic  equation  in  z.  Since  that  equa- 
tion must  hold  for  all  z the  coefficients  of  the  various  powers  of  z 
must  be  zero.  The  result  is 


C 3E 

r 


i 3_  (r  + 1_ 
r 3r  3r  ^ 2 


3^D, 


2 2 

30  36 


D2  (y+1)  = 0 


(66) 


(67) 


1 ^^^2 
J £ = n 

2 2^ 
30 


(68) 
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Also,  substituting  equations  (63)  and  (65)  into  equation  (34) 


gi  ves 


3r  30"' 


(69) 


Substituting  the  expressions  for  Vr,  v^,  and  v^J  into  the  four 
boundary  conditions  will  yield  additional  relationships  between  the  E's 
and  D's.  Applying  boundary  condition  of  equation  (60) 


Cl  + E, 


0,0 


3E, 

+ — - 

dr 


z = 0 


'0,0 


(70) 


Since  this  expression  must  hold  for  all  values  of  0 and  z. 


^1 

= 0 

(71) 

'll 

1 = 0 
*0,0 

(72) 

9Eoi 

C \ 

Or  1 

o 

n 

o 

o 

(73) 

Applying  boundary  condition  of  equation  (61)  and  setting  the  coefficients 
of  powers  of  z to  zero  gives 


' 0=fi/2 


Finally,  applying  boundary  condition  of  equation  (59),  and  setting 
the  coefficients  of  the  powers  of  z to  zero  yields. 


h2  30 
0 


1 9E2 


, 3h  3D,| 

1 _o  __2  = n 

h^  ao  30  . 

0 'h^. 


The  next  step  in  deriving  a solution  is  to  choose  a form  for  E, 
and  E2.  Again,  the  series  technique  is  employed  and  general  rela- 
tionships depending  on  r and  9 are  assumed  which  will  result  in  a re- 
duction to  the  axisymmetric  case  when  the  theta  dependence  is  surpressed. 
Then,  let 

E2  = D3(0)r  + 04(e)  (83) 

where  0^  and  0^  are  functions  of  0 only.  But  to  satisfy  equation  (73), 
04(0)  = 0 (84) 

Thus  after  integration, 

£2=03^+  DgO)  (85) 

Equations  (63),  (64),  and  (65)  indicate  that  D5  can  be  set  to  zero  since 
it  can  be  combined  with  0]  in  equations  (64)  and  (65)  or  is  zero  when 
di fferentiated  with  respect  to  r in  equation  (63).  Then, 


^2  ~ ^3  2 


Next,  integration  of  equation  (68)  twice  yields  D2  = C?-  + C3 
where  C2  and  C3  are  constants.  But,  to  satisfy  equation  (82)  C2  = 0 
and 

D2  = C3  (87) 

Substitution  of  the  above  expressions  for  and  D2  into  equation  (67) 
gives  the  following  partial  differential  equation: 

1 3 D3  2 1 ^ 

2D3  + ^ —2“  = (Y+DC3  - V (88) 

Using  the  technique  of  series  solution  once  more,  the  function  E-] 
is  assumed  to  be  of  the  form, 

El  = D7(0)r^  + Dg(9)r  + Dg(6)  (89) 

where  05,  D7,  Ds,  and  Dg  are  functions  of  0 only.  Boundary  equation 
(72)  indicates  that 

09(0)  = 0 (90) 

Substituting  equations  (88)  and  (89)  into  equation  (66)  gives 

2 Eg  1 3E 2 

r 40g  - ^Y+1)D3  + 3rD7  + (2Dg  - C3(y+1)D.|)  + ^ = 0 (91) 

— _J  y* 


In  addition,  applying  equation  (89)  to  equation  (69)  and  integrating, 
3Dg  .4  30.  3 3Dg  2 

3“  ^ 2"  ^10^®^ 


where  D]q  is  a function  of  0.  Differentiating  Equation  (92)  with  respect 
to  0 and  substituting  into  equation  (91)  yields 


2 2 

30 


1 

A] 

■4 

30  > 

1 

^*^10 

2 

30 

4Dj  - 2^r*l  )D3  t 1 t r[3D,  + 1 - Jj,  [2O3  - 


which  must  hold  for  all  values  of  r.  Thus 
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(94) 


(95) 

(96) 

(97) 

(98) 


where  is  a constant. 

Inspection  of  the  above  group  of  equations  shows  that  in  order  to 
solve  differential  equations  (94)  and  (96),  forms  for  D3  and  0]  must  be 
assumed.  Since  Dq  is  a function  only  of  0,  and,  since  the  flow  in  the 
tangential  direction  is  periodic,  the  function  D3  can  be  represented  as 
a Fourier  series  ( ^0) 


where  Cr,  Lp,  Mp,  and  Up  are  constants.  Furthermore,  in  order  that 
equation  (88)  be  satisfied,  Di  must  also  have  a Fourier  series  repre- 
sentation. 

00  ^ V OO  f ' 

“1  • Ce  * * ^1 

n=  I n=  I ' 

where  C5,  N^,  Rp,  and  6^  are  constants. 


^^^Spiegel,  M.R.,  Fourier  Analysis,  Schaum's  Outline  Series,  McGraw- 
Hill  Book  Co.,  1974. 


Since  can  be  written  in  terms  of  by  equation  (86),  bound- 

dO  f3U 

ary  conditions  (74)  and  (77)  can  be  used  in  evaluating  some  of  the 
constants  in  equation  (99).  That  is. 


dD, 

aF 


'0=0 


y u„L„cos 
n n 
n=1 


anO 


) M^(x„sin  « e 

n^l  n n In 


'0=0 


(101) 


The  re fore , 

Also, 

3D, 

(hT 


Ln  = 0 


(102) 


. 

00  , 

' 6=77/2 

n=  1 

1 J 

j 

= 0 


6=tt/2 


(103) 


Therefore, 


“n  = 


(104) 


Similarly,  boundary  conditions  (74)  and  (77)  applied  to  equation 
(100)  yields: 


I 

i 


Nn  = 0 

6n  = 2n 

Thus,  and  D^  are  written  as: 

oo 

D,  = y R cos2n0  + 

1 n=l  " ^ 

00 

0,  = y M„cos2n9  + Cr 
3 n 5 


(105) 

(106) 


(107) 


(108) 


The  above  two  expressions  may  now  be  substituted  into  equation 
(88)  to  obtain: 


y M„cos2n0 
n=l  " 


2 - iiDi: 
2 


1 * ? 9 

- ^ I 4n'^R^cos2ne  + 26^  - (y+1)C3=  0 (109) 

r n=l 
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Since  this  expression  must  hold  for  all  values  of  r and  0,  the  second 
two  terms  give  the  following  information: 


R = 0 
n 


^5  2 ^3 


In  addition,  it  is  seen  that  in  order  to  have  the  first  term  equal  zero 
for  all  0,  n = 1.  Then,  equations  (107)  and  (108)  reduce  to 


^1  = Cg 


= M^cos20  + C3 


Expressions  for  both  D]  and  D3  are  now  known  so  that  equations 
(94)  and  (96)  can  be  solved.  Equation  (94)  is  written  as 


16Df,  = 2C^(Y+l)fMTCOs20  + C3I 


Then,  a particular  solution^^^^  is  given 


(e-y)rp  (e-y)r, 
e -e 


g(y)dy 


where  r-j  = 4i , r2  = -4i , and 


g(y)  - 2C3(y+1  )[MiC0s2y  + C3 


Carrying  out  the  indicated  integration  and  adding  the  complementary 
solution  to  D5P  gives  the  complete  form  for  D5: 


Bellman,  R.,  and  Cookt , K.L.,  Modern  Elementary  Differential 
Equations , Addison-Wesley  Publishing  Co.,  Inc.,  1968. 


Dg  = CySin40  + C^cosAO  - - - ^ (cos40  - cos2o) 


"'16"-  - 1) 


where  and  Cg  are  constants. 

Looking  to  equation  (95),  its  solution  is  simply 


= CgSin3o  + CigCOs3'l 


where  Cg  and  C-|q  are  constants. 

Similarly,  the  solution  to  equation  (96)  is  seen  to  be 


Dg  = C-|^sin20  + C^2COs2o  + 


where  Ci],  C12.  and  C,g  are  constants.  Upon  substituting  this  expres- 
sion for  Dg  back  into  equation  (96),  C^g  is  found  to  be 


Then , 


^13  - ' 2 


Dg  = C^^sin20  + C^2COs20  + -y—  (r*!) 


Expressions  for  D5,  D7,  Dg,  and  Diq  are  now  known  in  terms  of  0 
and  various  constants.  Some  of  these  constants  can  be  evaluated  from 
boundary  conditions  (76)  and  (79).  First,  equations  (98),  (117),  (118), 
and  (121)  are  substituted  into  equation  (92),  the  result  being 


J r 2(Y+l)CgM,  (Y+1) 

Eg  = y 4CyCOs40  - 4CgSin40 -^2 (-4sin40  + 2sin20)  - — ^ 


"1  r^  r^ 

(-4sin40)  [+  j-  (3CgCOs30  - 3C^gSin30)  + ^ (2C^^cos20  - 
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When  condition  (7b)  is  applied  above. 


Cy  - Cg  - C-j^  - - 0 (123) 

Likewise,  applying  condition  (79),  it  is  found  that 


O 

o 

II 

o 

(124) 

Thus  far  the  boundary  equations  (80)  and  (81)  have  not  been 
satisfied,  A summary  of  the  results  thus  far  follows: 

v^  = El  + zrDg 

(125) 

= Cj  * Cf. , o/n 

(126) 

Vq  = Eg/r  + (DD2/30)(zr/2) 

(127) 

where 

t. 

El  = D^r^  + Dgr 

(128) 

Eg  = (3Dg/30)(r'^/4)  + ( 3Dg/36)  ( r^/2) 

(129) 

1 

Dg  = MiC0s26  + Cg  ((y+1)/2) 

(130) 

Y+1  (y+i)^Cg 

Dg  = Cgcos40  - C3M1  (cos40-cos20)  - — -^jg— ^(cos40-l ) 

(131) 

1 

and 

1 

^8  ^ C12COS20  + ^ CgCg 

(132) 

1 

All  that  remains  is  to  determine  the  constants  Cj,  Cg,  Cg,  C12 
and  M]  from  the  remaining  boundary  conditions.  The  nrst°boundary 
condition  to  be  used  is  equation  (81).  From  Equation  (41), 

1 

= h(O,0)  = Bq/[b^cos^0  + A^sin^0 

1/2 

(133) 

3h„  -h^  [a^  - bJI  sin20 

0 _ 0 t 0 oj 

30 

(134) 

i 
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Also,  from  equations  (43)  and  (58), 


m 


I I 


3 2 

cos  0 


kA^^sin^G 


(135) 


After  equations  (133),  (134),  (135),  (112),  and  (86)  are  substituted 
into  equation  (81)  and  evaluated  at  ho,  all  resulting  trigonmetric 
functions  can  be  written  in  terms  of  even  powers  of  cos6.  Each  re- 
spective power  of  cosO  is  factored,  leaving 


2 

cos  0 


^ ^ - A^l 


9 kA„ 


= 0,  (136) 


and 


C3  A ^ - A ^ = 0. 

3 2 0 0 1 B 

0 


(137) 


Since  these  two  expressions  must  be  true  for  all  values  of  0,  equation 
(136)  yields 

"1  - W <®o  - “''o> 


and  equation  (137)  yields 


B-W  («o  ^ 


1/2 


(139) 


I 


Following  the  same  procedure  with  equation  (80)  and  after  con- 
siderable algebraic  nenipulation  one  obtains 


. Ih  ( 


T6^-‘'T5^)-‘'T2T5T^ 


4T, 


^6  = 


2To(2T^5  " ^6) 


+ 


^^1^6  " ^^2^5  * ^^5^6^12 


(140) 

(141) 


h 


4 
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(142) 


where 


'oi6_  % 


2 

2 


Yiic  2 
2 ^3 

(143) 

(T^M^)/(3C3) 

(144) 

(1/7(403) 

(145) 

(146) 

(a/  - b/)/2b/ 

(147) 

In  summary  the  velocity  components  Vf.^,  v^3]  and  v^i  at  any  point 
0,  and  z can  be  easily  calculated  in  terms  of  the  geometric  variables 
Aq,  Bq.  and  k by  calculating  the  ten  constants  M-] , C3,  Tq,  Ti,  T2, 

^5’  ^6’  ^12’  ^6’  ^8  that  order  from  equations  (138),  (1  39), 

(143),  (i44),  (145),  (146),  (147),  (140),  (141)  and  (142)  respectively 
The  velocity  components  are: 

v^^=  cos40(Cg-T^-T2)r^  + cos20(Tir^  + C-|2r  + M^rz)  + 

.3  . ..  048) 

(^2'-  " " ^0^^) 


Vj^^=  -sin49(Cg-T^-T2)r^  - sin20(T^r^/2  + C^2^  049) 

v^^=  cos20(M^r^/2)  + T^r^/2  + Cgi  + Cg  (150) 


where  the  original  bar  notation  and  subscript  one  have  been  revived  to 
stress  that  the  result  is  nondimensionalized  and  is  only  the  first 
order  solution. 
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SECTION  V 


REDUCTION  TO  THE  AXISYMMETRIC  CASE 


In  the  first  appendix  of  reference  J the  analogous  solution  for 

the  flow  field  in  a converging-diverging  axisymmetric  nozzle  is  derived. 

A necessary  condition  for  the  elliptic  solution  of  this  investigation  is 

that  it  reduces  to  the  axisymmetric  solution  when  A = B and  k = 1. 

0 0 

The  first  order  solution  of  reference  3 for  non-dimensional  case 
(i.e.  Aq  = = 1)  is: 

u^  = y^/2  - 1/4  + [2/(y+l)]^/22 

(151) 

V-,  = (y/4)(y^-l)[(Y-H)/2]^/^  + yz 

(152) 

For  axisymmetric  flow  (i.e.  A = B and  k = 1)  equations  (138) 
through  (147)  reduce  to:  ° 

^1  ^1  " ^6  ' ^12  ■ ° 

(153) 

C3  = [2A^/(y  + D]^^^ 

(154) 

T2  = a//(4C3)-  Tj.I 

(155) 

Ce  ■ -V‘''  = '2  = 

(156) 

Thus 

for  axisymmetric  flow  equations  (148)-(150)  become 

Vi  = (r/4)(r^-l)[(Y  + l)/2Aj^/2A/  + A^fi 

(157) 

^01  = 0 

(158) 

^z1  = - Aq/4  + [2Aq/(y+1)]^/2 

(159) 

Equations  (157)-059)  are  identical  to  equations  (151)  and  (152)  (except 
for  differences  in  notation)  when  A^  is  set  to  unity. 
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SECTION  VI 


SUMMARY  AND  RECOMMENDATIONS 


1.  The  first  order  approximation  to  the  transonic  flow  field  for  a con- 
verging-diverging propulsion  nozzle  having  an  elliptical  cross  sec- 
tion has  been  obtained. 

2.  The  method  presented  is  generalized  such  that  it  can  be  extended  in 
a straight  forward  manner  to  obtain  higher  order  approximations  to 
the  flow  field. 

3.  The  solution  is  periodic  in  the  tangential  velocity  component,  vq, 
as  would  be  expected,  and  the  solution  reduces  to  the  special  case 

of  axisymmetric  flow  when  terms  containing  V0  and  8/36  are  surpressed. 

4.  The  solution  can  be  easily  programmed.  Parametric  studies  of  the 
flow  field  as  a function  of  geometric  variables  such  as  the  ratio  of 
major  to  minor  axes  of  the  elliptic  cross-section  could  provide  use- 
ful information  for  advanced  propulsion  nozzle  design. 
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